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Abstract 

A quantum algorithm for approximating efficiently 3-manifold topological in- 
variants in the framework of SU(2) Chern-Simons-Witten (CSW) topological 
quantum field theory at finite values of the coupling constant k is provided. The 
model of computation adopted is the g-deformed spin network model viewed as a 
quantum recognizer in the sense of 1 1 1, where each basic unitary transition function 
can be efficiently processed by a standard quantum circuit. 

This achievement is an extension of the algorithm for approximating poly- 
nomial invariants of colored oriented links found in [2 3 |. Thus all the signif- 
icant quantities -partition functions and observables- of quantum CSW theory 
can be processed efficiently on a quantum computer, reflecting the intrinsic, field- 
theoretic solvability of such theory at finite k. 

The paper is supplemented by a critical overview of the basic conceptual tools 
underlying the construction of quantum invariants of links and 3-manifolds and 
connections with algorithmic questions that arise in geometry and quantum gravity 
models are discussed. 
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1 Introduction 



The possibility of computing quantities of topological or geometric nature was recog- 
nized as a major achievement for quantum information theory in a series of paper by 
Michael Freedman and co-workers |4][5]|6) ( see a ^ so G) f° r an introduction). Their 
'topological quantum computation' setting, intrinsically fault-tolerant and protected 
from decoherence, was designed to comply with the behavior of 'modular functors' 
of 3D Chern-Simons-Witten (CSW) non-abelian topological quantum field theory 
(TQFT) [8] [9] ED* the gauge group being typically SU(2). In physicists' language, 
such functors are partition functions and correlators of the quantum theory and, ow- 
ing to gauge invariance and invariance under diffeomorphisms, which freeze out local 
degrees of freedom, they share a global, 'topological' character. More precisely, the 
physical observables are associated with topological invariants of knots -the prototype 
of which is the Jones polynomial ifTTl - and the generating functional is an invariant of 
the 3-dimensional ambient manifold, the Reshetikhin-Turaev-Witten invariant |[T2l l8l. 

The search for efficient quantum algorithms to compute (approximations of) knot 
invariants has been carried out by several groups in the last few years. 
Within the framework of topological quantum computation the existence of an effec- 
tive procedure has been taken for granted by resorting to the fact that this model is 
polynomially reducible to the standard quantum circuit model lfl3l . Later on, this im- 
plicit proof has been supported by the fundamental notion of 'additive approximation' 
introduced in lfT4l and borrowed in all the alternative approaches. 
The first explicit algorithm for the Jones polynomial of a knot presented as the plat 
closure of a braid has been constructed in lfl5l by efficiently approximating unitary 
matrices associated with a particular representation of the braid group, while in |[T6l 
the knot presentation and the representation of the braid group are generalized (these 
notions will be defined in the following sections). It is worth stressing that both these 
approaches rely on the standard model of quantum computation, namely the quantum 
circuit model based on qubits, elementary unitary gates and related algorithmic tech- 
niques such as the so-called Hadamard trick ifTJl . 

The approach we proposed in [2, 3| differs from the previous ones in many respects. 
First, we were able to handle more general knot invariants, namely 'colored' link poly- 
nomials (a link is a multicomponent knot) which reduces to Jones' in a particular 
case. Secondly, the model of quantum computation we adopted, that we refer to as 
the g-deformed spin network model, is in some sense a blending of the computational 
schemes mentioned above framed within the background provided by the theory of 
finite-states quantum machines This topic will be addressed in Section 2, while 
a few conceptual questions and implications of our model of computation will be dis- 
cussed in the last section of this introduction. 

In the main part of the paper, section 3, we extend our efficient algorithm for ap- 
proximating colored link polynomial to deal with 'quantum' 3-manifold invariants. 

In order to overcome the difficulties due to the fact that we shall resort to concepts 
and definitions arising in many different contexts (low-dimensional topology, quantum 
group theory, CSW field theory, 3-dimensional quantum gravity models, classical and 
quantum complexity theory), we shall illustrate in the next few paragraphs the basic 
conceptual tools underlying the construction of such 'universal' invariants. 
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1.1 From quantum topology to topological quantum field theory 



The term 'quantum topology' was introduced by Turaev [ 17 1 to denote implications on 
the topological side of the algebraic theory of quantum groups -technically, deforma- 
tions of universal enveloping algebras of Lie groups. The latter, based on the pioneering 
work of Drinfel'd 1 18 1 and Jimbo [ 19 1, was inspired by theoretical physics from its very 
beginning since quantum groups and associated i?-matrix representations are the basic 
tools of quantum inverse scattering methods and are the backbone of exactly solvable 
models in statistical mechanics 11201 . 

The deformation parameter q was originally assumed to be a real number related 
to Planck constant by q = e , therefore it is commonly referred to as a 'quantum' 
deformation, while the 'classical', undeformed Lie group simmetry is recovered at the 
particular value q = 1 (h — ► 0). When dealing with quantum invariants of knots 
and 3-manifolds ifTTl |2D q is most often a complex root of unity, the case q = 1 
being considered as the 'trivial' one. However, in a topologist's language, 'classical' 
topological invariants are not the (q = l)-counterparts of 'quantum' invariants, but 
rather the usual invariants of algebraic and geometric topology, typically related to the 
fundamental group and homology groups of manifolds and submanifolds. 

As Roberts remarks in the introduction to ||2T1 the standard topological invariants 
were created in order to distinguish between things and, owing to their intrinsic defi- 
nitions, it is clear what kind of properties they reflect. For instance, the Euler number 
X of a smooth, closed and oriented surface S determines completely its topological 
type and can be defined as x(S) — 2 — 2g, where g is the number of handles of S. 
On the other hand, quantum invariants of knots and 3-manifolds were discovered, but 
their indirect constuction based on quantum group technology often hides information 
about the purely topological properties they are able to detect. 

What is lost at the topological level is however well paid back by the possibility 
of bridging this theory with a plenty of issues in pure mathematics and theoretical 
physics (cfr. the review [21 1 and the list of references therein). To the early connections 
mentioned above (quantum inverse scattering, exact solvable models) it is worth adding 
the operator algebra approach used originally by Jones in defining his knot polynomial 
IfTTl . However, the most profitable development of the theory was that suggested by 
Schwarz and formalized by Witten [8 1 (see ll22l for a review and original references). 

Indeed, recognizing quantum invariants as partition functions and vacuum expec- 
tation values of physical observables in Chern-Simons- Witten topological quantum 
field theory provides a 'physical' explanation of their existence and properties. Even 
more radically, one could speak of a 'conceptual' explanation, as far as the topologi- 
cal origin of these invariants keeps on being unknown. In this wider sense, quantum 
topology might be thought of as the mathematical substratum of an SU (2) CSW topo- 
logical field theory quantized according to the path integral prescription (the coupling 
constant k > 1 is constrained to be an integer related to the deformation parameter q 
by q = exp(^|)). 

The CSW environment provides not only the physical interpretation of quantum in- 
variants but it does include as well all the historically distinct definitions 11231 . In 
particular, monodromy representations of the braid group [24| appear in a variety of 
conformal field theories since point-like 'particles' confined in 2-dimensional regions 
evolve along braided worldlines ( l25ll and references therein). As a matter of fact, 
the natural extension of CSW theory to a 3-manifold Ai 3 endowed with a non empty 
2-dimensional boundary dAi 3 induces on dAi 3 a specific quantized boundary con- 
formal field theory, namely the SU (2) Wess-Zumino- Witten (WZW) theory at level 



3 



I — k + 2 [8, 26 1 . The latter provides in turn the framework for dealing with SU(2) q — 
colored links presented as closures of oriented braids and associated with Kaul unitary 
representation of the braid group ll28l l27l . A further extension of this representation 
proposed by the same author in ll29l is used in this paper to construct explicitly the 
quantum 3-manifold invariants in the form originally defined in [ 30 1 within a purely 
algebraic setting. Such quantities are essentially the Reshetikhin-Turaev-Witten in- 
variants [ 12 1 evaluated for 3-manifolds presented as complements of knots/links in the 
3-sphere S 3 , up to an overall normalization. 



1.2 Algorithmic complexity of Chern-Simons-Witten theory 

As mentioned in the introductory remarks, the 'quantum field' computer |4l|5]|6) is a 
model of computation designed to process quantities of topological nature arising in 
CSW environment and thus 'effienciency' of any calulation -such as that of Jones knot 
polynomial- should be guaranteed by definition (we leave aside here the issue of 'ex- 
act' versus 'approximate' calculation at least for the moment). 

However, when dealing with algorithmic questions, the model of computation adopted 
should comply with the commonly accepted paradigms of theoretical computer sci- 
ence. Turing machines, together with the polynomially equivalent circuit models based 
on elementary boolean gates, represent the universal schemes which allow problems 
and algorithms to be grouped into classical complexity classes BTI . In quantum com- 
puting, the notions of quantum Turing machine -and associated quantum circuits based 
on qubits and unitary elementary gates- can be introduced and represent the standard, 
universal model of computation lfl3l . 

It was shown in [5] that SU (2) CSW functors at the fifth root of unity, whose domain 
is restricted to collections of 'topological' qubits (disks with three marked points) on 
which suitable unitary representations of the braid groups B3 and B3 x B3 act, do 
reproduce the standard elementary gates of the quantum circuit model. The physical 
support of such information processing consists of anyonic systems obeying particular 
types of braid statistics, and work is in progress on the experimental side to check the 
implementability of such approach (see |32| and also [33] for a different kind of im- 
plementation of braiding operators). 

Based the above properties, a sort of 'minimal' realization of the full quantum field 
computer is (polynomially equivalent to) the quantum circuit, and indeed the work of 
ifTSl and lfl6l on quantum computation of the Jones invariant does not depend at all on 
any field theoretic background. Generally speaking, this is satisfactory for the aim of 
introducing an ad hoc computing model for treating anyonic quantum systems; yet is 
somehow disturbing because classical Turing machines and their probabilistic counter- 
parts are able to simulate efficiently only any dynamical systems governed by classical 
laws at any degree of accuracy. Even the objection that we are in the presence of a 
quantum field theory -not simply a quantum mechanical many-body system- is mis- 
leading since the CSW model is exacly solvable at the full quantum level (for each 
fixed value of the coupling constant) without resorting to any approximation such as 
the low-energy limit l22l [231 . The crucial feature of possessing only global, purely 
topological degrees of freedom makes quantum CSW theory likely to be simulated 
within a computational scheme based on a discrete space of states and able to imple- 
ment polylocal braiding operations. 

As will be illustrated in section 2, the universal model of computation able to handle 
all discrete, many-body quantum systems described by (real or virtual) pure angular 
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momenta states (not simply two-level systems) as well as solvable field theories of 
CSW-type is the (q-deformed) spin network proposed in (2] [3]. 

As already recognized in [4 1 the task of computing the Jones polynomial of a knot 
or link represents a major achievement since it is the simplest observable in quantum 
CSW theory and then its calculation is a testing ground of the effectiveness of topo- 
logical quantum computation. However, the interest in this algorithmic problem has 
recently grown in connection with the search for new testing grounds for quantum in- 
formation theory in general, without necessarily exploiting the physical meaning of the 
invariant. 

The reason why Jones link polynomial is so crucial in the computational context 
relies on the fact that a 'simpler' link invariant, the Alexander-Conway polynomial, 
can be computed efficiently, while the problem of computing 2-variable polynomials 
-such as the HOMFLY invariant- is NP-hard (see [34] [21] [35] for definitions of these 
invariants, original references and accounts of algorithmic questions). 

The issue of computational complexity of the 1-variable Jones polynomial in clas- 
sical information theory can be summarized in 

Problem 1 

How hard is it to determine the Jones polynomial J{L, q) of a link L? 

A quite exhaustive answer has been provided in (JVW), where the evaluation of the 
Jones polynomial of an alternating link L at a root of unity q is shown to be ^P-hard, 
namely computationally intractable in a very strong sense. 

Recall first that 'alternating' links are special instances of links, the planar diagrams 
of which exhibit over and under crossings, alternatively. Thus the evaluation of the 
invariant of generic, not only alternating links is at least as hard. Secondly, the compu- 
tation becomes feasible when the argument q of the polynomial is a 2nd, 3rd, 4th, 6th 
root of unity, so that the first not easy case involves a 5th root of unity (refer to ll35l 
for details on this technical issue). Finally, the jfP complexity class can be defined 
as the class of enumeration problems in which the structures that must be counted are 
recognizable in polynomial time. A problem tt in #P is said to be #P-complete if 
for any other problem it 1 in #P, tt' is polynomial-time reducible to ir; if a polynomial 
time algorithm were found for any such problem, it would follow that #P C P. A 
problem is #P-hard if some #P-complete problem is polynomial-time reducible to 
it. Other instances of #P-complete problems are the counting of Hamiltonian paths 
in a graph ||3T1 and the more intractable problems arising in statistical mechanics, such 
as the enumeration of all configurations contributing to ground state partition functions 

The computational intractability of Problem 1 does not rules out the possibility of 
'approximating' efficiently Jones invariant, so that we may ask 

Problem 2 

How hard is it to approximate the Jones polynomial J(L, q) of a link L at 
a fixed root of unity q (q ^ 2nd, 3rd, 4th, 6th root)? 

Loosely speaking, the approximation we are speaking about is a number Z such that, 
for any choice of a small r\ > 0, the numerical value of J(L,q), when we substitue 
in its expression the given value of q, differs from Z by an amount ranging between 
— r\ and +77 (see section 3.2 below for a more precise statement). In the framework of 
classical complexity theory no algorithm to handle efficiently Problem 2 exists, while 
the answer in the quantum computational context was given in 1 14| (see also 1 16|): 
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The approximation of the Jones polynomial of a link presented as the clo- 
sure of a braid at any fixed root of unity is BPQ-complete. Moreover, this 
problem is universal for quantum computation, namely is the 'prototype' 
of all problems efficiently solvable on a quantum computer. 

Recall that BQP is the computational complexity class of problems which can be 
solved in polynomial time by a quantum computer with a probability of success at 
least | for some fixed (bounded) error. In lfT4l it was proved that P J = BQP, where 
P J is defined as the class of languages accepted in polynomial time by a quantum 
Turing machine with an oracle for the language defined by Problem 2. This equality 
between computational classes implies that, if we find out an efficient quantum algo- 
rithm for Problem 2, then the problem itself is complete for the class BQP, namely 
each problem in this class can be efficiently reduced to a proper approximate evaluation 
of the Jones polynomial of a link lfl6l . 

Explicit, efficient quantum algorithms for approximating the Jones polynomial 
were proposed by Aharonov, Jones and Landau [15] and by Wocjan and Yard lfl6l . 
while an early attempt can be found in |36|. 

In our papers flU O we proved that efficient algorithms can be implemented for ap- 
proximating the larger class of 'colored' Jones polynomials of links (addressed also in 

E3). 

The issue of colored link invariants brings us back to the quantum CSW environ- 
ment, where they represent the most general gauge invariant physical observables of the 
theory, being vacuum expectation values of generic Wilson loop operators l8l l22l l23ll . 
In section 3 we shall provide a generalization of our algorithm for colored polynomials 
to handle the quantum 3-manifold invariants introduced in 1.1. 

Summarizing our results, we have shown that all the significant quantities -partition 
functions and observables- in SU (2) quantum CSW theory can be efficiently approx- 
imated at finite values of the coupling constant k. The intrinsic field-theoretic solv- 
ability of CSW theory is thus reflected by its computability on a quantum computer. 
Looking at the question the other way around, classical computational intractability 
of Jones and colored polynomials can be viewed as a consequence of their quantum 
nature. This feature has prevented up to now both exact and approximate efficient 
calculations of such topological quantities on classical (probabilistic) machines as it 
happens for simulations of any non trivial 'genuine' quantum mechanical system (see 
e.g. Feynman's proof in |38|). 

1.3 Quantized geometry versus quantum computing 

In this section we address some implications of our results in connection with algorith- 
mic questions that may arise in other physical theories whose dynamical variables have 
a geometric character, typically quantum gravity models in D = 3 and 4 spacetime di- 
mensions. 

We have however to face preliminarly a conceptual dilemma, namely whether 

i) an abstract universal model of computation, able to simulate any discrete 
quantum system including solvable topological field theories, must exists 
by its own 

or 
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ii) a (suitably chosen) quantum system is by itself a computing machine 
whose internal evolution can reproduce the proper dynamics of classes of 
physical systems. 

The second alternative is becoming quite popular thanks to Lloyd's model, where a net 
of computing units generates a (superposition of) AD spacetimes ll39l (see also ll40l 
where similar ideas were anticipated). 

The idea that many (if not all) aspects of our reality may be thought of as 'outputs' 
of some kind of information processing is both appealing and intriguing. In this con- 
nection the role of information theory and its tools is so enhanced that it becomes a 
unifying paradigm. Of course the classical version of hypothesis i) is usually taken for 
granted as far as, on the one hand, a (probabilistic) Turing machine is capable of simu- 
lating the evolution of any classical system within a given accuracy, and, on the other, 
all concrete, finite-size realizations of the abstract machine obey the laws of classical 
physics. The praticability of hypothesis ii) depends heavily on which system is chosen 
as a simulator and which types of boundary or initial conditions must be imposed to 
reproduce the dynamical behavior of observed physical systems. Moreover, the con- 
cept of 'efficient' processing of information seems difficult to be handled without an 
abstract reference model of computation. 

With the previous remarks in mind, we favour assumption i), where the reference 
model of computation can be represented by the spin network simulator ATI l42l [431 
(see also [44]) or its q-deformed, finite-size counterpart ||2][3) (see also section 2 be- 
low). Note however that it may be tempting to assert that the spin network -thought 
of as a real net of interacting spin variables- can play as well the role of the reference 
quantum system in statement ii). This is due to the fact that the recoupling theory of 
SU(2) angular momenta -representing the algebraic substratum of the simulator- is 
the main ingredient of the 'spin network models' introduced by Ponzano and Regge 
ll45l and Penrose [46 1 . Here classical, discretized 3D euclidean geometry arises from 
recouplings of quantum angular momenta in the asymptotic, large angular momentum 
limit (see ll42l . section 5 for a brief account). Thus, much in the sense of ii), spin 
networks may act -under suitable constraints- as computing machines able to process 
information encoded into quantum spins to output 'quantized' 3-geometries obeying 
Einstein equations in the classical limit. 

Thinking back to the issue of algorithms for quantum invariants defined in the 
framework of quantum CSW theory, our results can be used to test the algorithmic 
complexity of quantum gravity models in D = 3 spacetime dimensions too. This 
achievement can be justified -independently from spin network models- by exploiting 
the close connection between CSW theory and 3D gravity both as classical field theo- 
ries and at the quantum level j8l|26] (classical 3D gravity with a positive cosmological 
constant, reinterpreted as an SU(2) gauge theory, is quantized through the Euclidean 
path integral prescription). 

For a closed orientable Riemannian 3-manifold M 3 , let Z(A4 3 ; k) denote the Witten 
partition function associated with the classical CS action Scs(A), A being the con- 
nection 1-form. The functional Z(M 3 ; k)Z(M 3 ; k) = \Z(M 3 ;k)\ 2 = Z(M 3 ;k) 
for finite k is the partition function of 3D Euclidean quantum gravity in the first- 
order form, where the coupling constant k is related to the cosmological constant A by 
£; = 47T/VA.Q 

'Note in passing that the invariant Z(M 3 ; k), on the one hand, equals the Turaev-Viro invariant for tri- 
angulated 3-manifolds (47'] and, on the other, can be derived by relating 3D gravity to a BF-type topological 
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The point here is that every manifold in the class considered here can be presented 
as the the complement of a link L in the 3-sphere (technically, by surgery along a 
framed version of L, see section 3.1), so that A4 3 = M\ ~ (S 3 \ L) and Z(M,\\ k) 
= \Z{M.\] k)\ 2 , where Z(M.\\ k) is the Reshetikhin-Turaev-Witten quantum invari- 
ant introduced in 1.1. Thus the results on algorithms for quantum invariants of links 
and 3-manifolds summarized in 1.2 work equally well in the context of Euclidean 3D 
quantum gravity models, reflecting once more the exact solvability of the theory for 
finite k. 

An even more interesting connection between SU(2) q quantum invariants and 
quantum gravity emerges when dealing with canonical quantization methods applied 
to general relativity in (3 + 1) dimensions. We refer in particular to the 'loop rep- 
resentation' based on Ashtekar 'connection representation' of canonical gravity (see 
Il48ll49 l for reviews and original references). The Reshetikhin-Turaev-Viro invariants 
act there as quantum states associated with the boundary 3-geometries (spatial slices 
of 4D spacetimes). Such states were shown to satisfy both the quantum diffeomor- 
phism constraint and the quantum Hamiltonian constraint (Wheeler-DeWitt equation) 
once a cut-off -related to the deformation parameter q- has been fixed. Thus quantum 
invariants constitute the natural kinematical arena of loop quantum gravity 11501 . This 
sketchy scenario hides however the necessity of introducing new types of quantum in- 
variants -the perturbative invariants- which arise when the cut-off of the underlying 
quantum field theory is removed. 

In the CSW field-theoretic setting, perturbative invariants emerge as coefficients 
of the asymptotic expansion of the partition function Z(M.\] k) as k — > oo (see ED . 
ch. 7). Besides a first term that corresponds to the semiclassical approximation of 
Z{M. L ]k) (saddle point), each contribution in the expansion is exp{2irkScs{A)} 
times a power series in 1/k (here Scs(A) is the CS action evaluated for the flat con- 
nection A). Perturbative invariants are the coefficients of the powers (1/fc)" evaluated 
by using (n + l)-loop Feynman diagrams. No complete perturbative treatment of 
quantum CSW theory is available at present, and the meaning of such invariants in ge- 
ometric topology is only conjectured in a few cases. 

We are interested here in discussing briefly the 'volume conjecture' concerning special 
classes of hyperbolic 3-manifolds (recall that the volume is a topological invariant for 
such manifolds). Note also that most manifolds obtained by surgery on framed links 
in the 3-sphere can be endowed with hyperbolic metrics. Focusing in particular on 
'hyperbolic knots', namely those knots which give rise to finite volume hyperbolic 3- 
manifolds, the volume conjecture proposed in BTl l52l (see also the review ||2T1 for 
extended versions) can be cast in the form 

^ ^ log\MK)\ = (1) 

where if is a hyperbolic knot and the notation Jjv (K) stands for the iV-colored poly- 
nomial of K evaluated at q = exp(2iri / N) . 

As pointed out many times, all quantum algorithms dealing with link polynomials are 
established for a fixed choice of the root of unity q appearing in the argument of the 
invariants, while the volume conjecture involves the analysis of the asymptotic behav- 
ior of single-colored polynomials of a same knot for increasing values of the coloring 

field theory |8 22 26]. 
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itself. 

It would be interesting to explore the possibility of borrowing some of the tech- 
niques employed in ll53l for dealing with the asymptotics of the Jones polynomial to 
test conjecture (HJ within the computational framework for colored polynomials ad- 
dressed in the present paper. 

Note finally that the relevance of this issue for quantum gravity stems from the ob- 
servation that all vacuum solutions of Einstein field equations in (2 + 1) dimensions 
with a negative cosmological constant are hyperbolic metrics. As recently shown by 
Carlip [54 1, the smallest hyperbolic volumes give rise to the largest contributions in 
the saddle point term of the path integral of the quantum theory. This implies in turn 
that the so-called 'real tunneling geometries' are most probable (a real tunneling ge- 
ometry represents the transition from a compact Riemannian spacetime to a Lorentian 
one within the framework of the Hartle-Hawking 'no boundary' approach to quantum 
cosmology [55 1). On the other hand, the possibility of controlling such invariants from 
the algorithmic point of view might help also in selecting weights to be assigned to 
3-geometries in the 'sum over topologies' within a 'fully quantum' cosmology theory. 

2 Processing braiding operators 

2.1 g-deformed spin network and quantum recognizers 

As pointed out in the introduction, our reference model of computation to deal with 
quantum topological invariants derived from Kaul unitary representations of colored 
oriented braids [27 , 28 1 is given by the q-deformed spin netwok model. 
The (undeformed) spin network simulator has been extensively addressed elsewhere 
[42 1 but for the convenience of the reader we have collected in appendix A a few 
mathematical details. However, in order to recognize the necessity of introducing a 
q-deformed version, it is worth to discuss here the (categorical) foundations of the 
quantum theory of angular momenta. The computational space of the spin network 
simulator -modelled as a graph- encodes the representation ring of the Lie group 
SU(2) -namely finite-dimensional Hilbert spaces supporting irreducible representa- 
tions (irreps) of SU(2) endowed with two binary operations, tensor product ® and 
direct sum © (that provide a ring structure over the field C)- together with all the 
unitary operators relating (multiple tensor products of) such spaces. Unlike the usual 
quantum circuit model ff3l . here it is possible to handle directly eigenstates of N 
binary coupled angular momentum variables labelled by integers and half-integers 
ji , j2 , . . . , j m (in h units) and not simply A-qubit states labelled by the fundamen- 
tal irrep ji = j2 = ■ ■ ■ = ]n = 1/2. The (re)coupling theory of N SU(2) angular 
momenta provides the whole class of unitary transformations that can be performed 
on many-body quantum systems described by pure angular momentum binary coupled 
sets of eigenstates ( ll56l . topic 12). The unitaries that we need here (referred to as 
j-gates in section 3 of ll42l ) are phase transformations, related to swaps of two con- 
tiguous spins, and recoupling trasformations expressed in terms of 3nj-coefficients of 
SU(2) and related to changes in the interaction schemes of the A angular momenta 
(A = n + 1). 

Within the framework of the categorical approach, the SU (2) representation ring 
is an instance of a unitary tensor category, endowed with intertwiner spaces and two 
basic morphisms, a 'twist' (a trivial type of 'braiding') and an 'associator' j57l . l58l . 
The former acts on the tensor product of two Hilbert spaces V, W supporting irreps of 
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SU(2) by exchanging the order of the factors, namely 

R v ,w ■ V <8> W -» W ® V (2) 

with Rw,v ° -Rv,w = Wy^vy- 

The explicit action of R on a (binary coupled) state is a trivial phase transform, see 
d47l >. (ijgT l in appendix A. 

The associator F relates different binary bracketing structures in the triple tensor prod- 
uct of irreps V, U, W 

F : (V <E> U) <E> W ^ V (g> (U (g> W) (3) 

and is implemented on a binary coupled state as a Racah transform involving one 
Racah-Wigner 6j-symbol (see d45l l. d46b in appendix A). Note that both and (01 
are isomorphisms (unitary morphisms between intertwiner spaces in the categorical 
language) but the associator reflects a true (physically measurable) modification of the 
way in which intertwiner spaces are coupled. 

The remarkable fact, derived from the general theory of braided tensor categories 
(59), is that more complicated multiple tensor products can be handled without intro- 
ducing any further independent morphism. Actually, multiple tensor product spaces 
can be related by different combinations of braidings and associators, so that the ba- 
sic morphisms must satisfy compatibility conditions, a so-called pentagon identity and 
two exagon identies. In the more concrete language of SU(2) recoupling theory, each 
3nj-symbol can be obtained as a combination of phase and Racah transforms and 
the non-uniqueness of such decomposition is translated into the Biedenharn-Elliott 
(pentagon) identity and the Racah identity (see e.g. Il56l or [60] for their explicit ex- 
pressions). 

In order to deal with non-trivial braiding operators -to be used in connection with 
the study of braid group representations and braid statistics (typically occuring in any- 
onic systems)- we are forced to modify the operation R in © by introducing a genuine, 
non-trivial braiding morphism |59| 

Tlv,w ■ V <g> W -» W ® V (4) 

with TZwy ° T^v,w Idy®w- 

A consistent way of modifying the SU (2) tensor category to include non-trivial braid- 
ings can be achieved by moving to the representation ring $H (SU (2) q ) of the q— 
deformed Hopf algebra of the Lie group SU(2), SU(2) q (q a root of unity). The 
resulting braided tensor category is the 'universal' algebraic structure underlying the 
constructions of quantum invariants of links and 3-manifolds outlined in 1.1 and 1.2. 
(We refer the reader to [2J and to [61 1 for short technical surveys of the quantum group 
and CSW approaches.). 

According to the above remarks, it should be clear that the 'g-deformed' spin net- 
work model of computation J2j[3| is modelled on the g-tensor category 

(m(SU(2) q ); TZ; T), (5) 

where we have denoted by T the (/-counterpart of the associator F in (O. Once suitable 
basis sets are chosen in the finite collection of irreducible spaces S s )\(SU(2) q ), the 
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unitary morphisms 1Z and T can be made explicit. In particular T turns out to contain 
the c^-deformed counterpart of the 6j-symbol and, regarding it as a unitary matrix, it is 
also referred to as 'duality' (or 'fusion') matrix borrowing the language of conformal 
field theories ll25l . 

The efficient quantum algorithm for (approximating) SU(2) q -colored link poly- 
nomials we obtained in [3] relies on a two-level procedure which can be summarized 
as follows. 

• Kaul unitary representation of colored oriented braids -associated with links pre- 
sented as plat closures of such braids- is processed on the q-spin network (0 in 
a number of steps that grows polynomially in the size of the input. In particular, 
each elementary computational step is implemented by applying either 1Z or T 
(see section 2.2 below). 

• The basic g-morphisms 1Z, T in the Kaul representation can be efficiently com- 
piled on a standard quantum computer, by means of universal elementary gates 
acting on suitable qubit-registers.|3 

In order to analyze in more details the first topic above we need to frame the q- 
deformed spin network model within the theory of quantum automata and quantum 
languages. According to [1 1 a quantum recognizer is a particular type of finite-states 
quantum machine defined as a 5-tuple {Q, H, X, Y, T(Y\X)}, where 

1 . Q is a set of n basis states, the internal states; 

2. H is an n-dimensional Hilbert space and we shall denote by l^o) £ H a start 
state expressed in the given basis; 

3. X and Y — {accept, reject, e} are finite alphabets for input and output symbols 
respectively (e denotes the null symbol); 

4. T(y|X) is the subset of n x n transition matrices of the form {T(y\x) = 
U(x)P(y); x <E X, y <E Y}, where U{x) is a unitary matrix which determines 
the state vector evolution and P{y) is a projection operator associated with the 
output measurement on (suitable complete sets of observables associated with) 
the upgraded state vector. 

In this kind of machine the output alphabet is chosen in such a way that a word w 
written in the input alphabeth X must be either accepted or rejected, while for the null 
symbol the requirement is P(e) = I (the identity matrix). Thus the one-step transition 
matrices applied to the start state |^ ) can in principle assume the forms 

a) T(e\x) = U(x)P(e) = U{x)l Vx G X, 

2 This result is quite interesting by itself as recently pointed out in [62'j where relations between the basic 
morphisms and large Fourier transforms are addressed. 

The problem of finding out efficient algorithms to compute the fundamental functions of the quantum theory 
of SU (2) angular momenta -Clebsch-Gordan coefficients, 6jr'-symbols- has attracted much attention 1631 . 
However, to our knowledge, there are neither classical nor quantum algorithms avalaible for evaluating these 
functions for arbitrary values of their arguments. The crucial remark is that a g-6j'-symbol with arbitrary 
entries can be efficiently compiled and approximated on a quantum circuit independently of the input size 
of the algorithmic problem owing to the presence of the natural cut-off provided by k (see section 3.2, in 
particular footnote 8). Then we realize once more that the (/-symmetry of solvable topological field theories 
is indissolubly tied with the effective computability of these models. 
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b) T(accept|a;) = [/(x)P(accept) Vi€l with P(accept) = |*o)(*o|, 

c) T(reject|a;) = [/(x)P(reject) \/x £ X with P(reject) = I- |*o)(*o|, 

according to whether no measure is performed (case a)), or the output is 'accept'/'reject' , 
namely cases b) /c) respectively. 

The general axioms stated above can be suitably adapted to make this machine able 
to recognize a language C endowed with a word-probability distribution p(w) over the 
set of words {w} £ C. In particular, for any word w — X1X2 ■ ■ ■ Xi G C the recognizer 
one-step transition matrix elements are required to be of the form Tij(x s ) = Uij{x s ) 
on reading each individual symbol x s £ w, namely no measurement is performed at 
the intermediate steps (here i, j run from 1 to n, the dimension of the Hilbert space 7i). 
Each Uij{x s ) must satisfy the condition 

\Uij(x s )\ 2 > 0, (6) 

and the recognizer upgrades the (normalized) initial state to 

U(w)\* ) = U(xi) . . . U{xi) | * > • (7) 

Then the machine assigns to the word w the number 

p{w) = | (* | U(w) P(accept) U(w)\^ ) | with < p(w) < 1 , (8) 

which corresponds to the probability of accepting the word w as a whole. 

More generally, the machine accepts a word w according to an a priori probability 
distribution Pr(w) with a word-probability treshold Q < S < 1 if 

\T?r(w) - p(w) | < 5, Vwe£. (9) 

In what follows the accuracy S will be set to 0, so that the two probability distributions 
Pr and p coincide. 

2.2 The Kaul representation as a quantum language 

In this section we shall show that the g-deformed spin network computational scheme 
(01 embodies families of quantum finite-states machines (or quantum automata) {A q } 
- parametrized by the labels of N irreps of SU(2) q , q =root of unity- that recognize 
the language generated by the braid group according to a probability distribution given 
by the square of the modulus of the g-colored link polynomial. 

This construction complies essentially with what we have done in [2], but here we 
should stress the interpretation of the Kaul representation 11281 as a quantum language, 
on the one hand, and the role played by the probability distribution, on the other (leav- 
ing aside details on the field-theoretic background material). 

The basic ingredients of Kaul's construction are oriented geometric braids (see Fig. 
H]), the strands of which are endowed with 'colorings' given by SU (2) q irreps labels. 
An n-strand colored oriented braid is defined by two sets of n assignments ji — (ji , ) 
with (i — 1,2, ... ,n), corresponding to the spin ji labelling the strand and to the ori- 
entation £j of the strand, with = ±1 (for the strand going into or away, respectively, 
from a horizontal rod from which the braid issues). The first set of assignments is 
associated to the upper rod, the second to the lower rod (we use the convention that 
two braids are composed in the downward direction). The conjugate of ji is defined as 
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Figure 1: An oriented braid on four strands. 

j* = (ji, — £,). It follows that the assignments on the lower rod are just a permutation 
of the conjugates of the assignments on the upper rod. A colored and oriented braid 
can thus be represented by the symbol 

o{\ f f Y (10) 

V n h ■ ■ ■ hi J 

where lj — j n u\ for some i and j and a permutation tt of {1, 2, . . . n}. 

The composition of two colored oriented braids is well defined only if the orien- 
tations and the colors of the two braids match at the merging points, as shown in Fig. 
. The group of colored oriented braid is generated by the identities (one for each 

IMI 

b, 



1 4 j l i 







M 















Figure 2: The composition of two colored oriented braids. 

assignment of colors and orientations on a topologically unentangled braid) and by the 
braids of type <ji, as shown in Fig. 0, 

x- 

fl % fn-1 J*" h P f f|-,1 fn1 f™ 

I CT| 

Figure 3: The identity and the generator ai of the colored oriented braid group. 
The collection of {07} for I — 1, 2, ... n — 1 are the standard generators of the 
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(colored) braid group B n and satisfy the following defining relations 



ai<Ji + i(Ji 



C"i+l<TiCTj+i 

CTjCTi \i - j\ > 2. 



(ID 



The inverse of a generator 07, (cr/) _1 , corresponds to the under-crossing of the left 
strand ji in Fig. (right). 

In order to obtain a link (multicomponent knot) from a colored braid we need to 'close 
up' the braid. For our purposes we may consider only the plat closure (or platting) of a 
colored braid, defined for braids which possess an even number of strands and whose 
assignments match as follows (see also Fig. J4) ) 



Ji 3i 32 32 
li l* I2 ^2 



32n 32 
hn I 



2n 



(12) 



A A* A A* A A* 
ji ji J2 J2 J2„ j2„ 



A A* A A* A A* 
ll ll I2 b bn bn 



platting 



. . A*A_ A* 

^ A A 



2n-plat closure 



*VfA/f 



I h b b bn I 



A/a* 



Figure 4: The platting of a colored oriented braid on 2n strands. 

Since any (colored oriented) link can be obtained as the plat closure of a braid [64|, 
we do not lose generality for what concerns the class of links that can be handled. 

The further step consists in embedding the 2n-strand braid into a 3-sphere 5* 3 
with two three-balls removed, giving rise to a 3-manifold with two boundaries E 1 , E 2 
(topologically two 2-spheres S 2 with opposite orientations). The intersections ('punc- 
tures') of the braid ( TT~2T > with the boundaries inherit the colorings and orientations from 
the corresponding strands of the braid (to be associated with Wilson line operators in 
the ambient CSW topological field theory). Following [8], finite dimensional Hilbert 
spaces Ti l ® Ti 2 are associated with to the two boundaries E , E 2 , and the basis sets in 
these spaces are the so-called conformal blocks of the boundary Wess-Zumino-Witten 
conformal field theory at level I, with 2n external lines labelled by (different) irreps of 
SU(2) q (I is related to CSW coupling constant by I = k + 2, so that from now on 
we set q — exp{27rz/£}). Two particular types of conformal block bases are needed to 
deal with braids the plat closures of which will give rise to colored oriented links, and 
their combinatorial patterns are shown in Fig. 13 and Fig. 151 . 




J2n-1 J2n 




l"l Vl fn-3 

Figure 5: The conformal block of type {p; r} (odd). 



The (orthonormal) basis sets are constructed by taking particular binary coupling 
schemes of the 2n 'incoming' angular momentum variables 31,32, ■■ ■, 32n which must 
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Figure 6: The conformal block of type {q; s} (even). 



sum up to give a spin-0 total singlet stated The procedure can be carried out by paral- 
lelling the SU (2), undeformed case (see appendix A) but here the labels of the irreps 
ji's (integer and half-integers) are constrained to range from to 1/2 and the binary 
bracketings on tensor products are decomposed according to rules of SU(2) q repre- 
sentation theory (see e.g. (6T|, section 3). 

Looking at the combinatorial structure of Fig. 0, the most general odd-coupled 
basis is consistently labelled as 

|p;r) j , (13) 

where j stands for the ordered string ji,j2, ■ ■ ■ J2n, P = Po,Pi, - ■ ■ ,Pn-i and r = 
r%, r2, ■ ■ ■ , r n _3. In the even-coupled case depicted in Fig. (6), the states of the basis 
are denoted by 

(14) 



|q;s, 



where j is the same as before while q = qo, qi, . . . , q n -i and s = si, S2, ■ ■ ■ , s n -3- 

The basis vectors associated to the conformal blocks ([L3l and <fT~4t are related to 
each other by 

h h 



p; r ; 



E 

(q;s) 



.4 



(q;s) 
(p;r) 



.73 



34 



J2n-1 J2r, 



q; s ) 



(15) 



where the symbol ^ip.*) [ : : ] represents the unitary duality matrix (or g-deformed 3nj 
recoupling coefficient). As pointed out in section 2.1, it is a standard result that any 
such duality matrix can be decomposed into (sums of) products of 'basic' duality ma- 
trices or q-6j symbols, see section 3.2 below. 

A graphical representation of the decomposition dT~5b in the case of eight incoming spin 
variables is shown in Fig. JT) . Note that the graphical representation of the basic du- 
ality transformation (the matrix counterpart of the associator in the language of tensor 
categories) can be also drawn in the most familiar form shown in Fig. (8). 

As proved in ESI , the colored polynomial of a link L, presented as the plat closure 

3l 3i ■ ■ ■ Jn j*, 

if* 11* 



a of a colored braid a 



as defined in (fT2l . is given by 



J[L-y, q] = Tr jj [2* + 1], \0;Q\U 



j=i 



Ji 3i 

ii it 



j f* 



|0; 0)J , 
(16) 



3 Incidentally, such a choice complies also with the conditions necessary for the construction of error 
avoiding codes and implies as well the robustness of the scheme 1651 1661 . 
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q-6j 



q-6j 



|q-6j 

PY| tf |YYY| 



Figure 7: Example of decomposition of a duality transformation between the two ex- 
tremal conformal blocks in the case of eight incoming spins. 




j, j. 



Figure 8: The elementary duality transformation associated with a q-Qj symbol. 



where j = (j%,j2, ■ ■ ■ ,j2n), + l]g is the quantum dimension of the irrep labelled 
by ji and the quantum integer [x] q is defined as 



r?/2 



--9/2 



r-1/2 



-1/2- 



(17) 



Thus J [L; j; q] can be evaluated by taking the trace of the matrix elements of the 
composite braiding operator U[a] in the Kaul representation with respect to the odd- 
coupled basis, where all the intermediate quantum numbers are constrained to give sin- 
glet eigenstates (a similar result would hold true for the even-coupled basis). Moreover, 
U[a] can be decomposed into a finite sequence of unitary matrices C/[(72/+i] (diagonal 
matrices in the odd-coupled basis adopted in ( TToT l) and duality matrices of the type (T5[ 
to be applied whenever a switch to the even-coupled basis is needed, namely when an 
even t/[<72;] occurs in the decomposition (see l28l for the explicit expression of these 
matrices). 

The construction outlined above can be cast into an effective process of calculation 
by resorting the concept of quantum recognizer introduced in section 2.1. 
The A q recognizer is defined, for a fixed root of unity q = exp{27ri /£}, by the 5-tuple 

{C odd , H, B 2n , {accept, reject, e}, U(B 2n )}, where 

• C dd is the odd-coupled conformal block basis of the boundary WZW theory 
(see (O and Fig. 0). 

• H is the ordered tensor product of 2n (2ji + 1 )-dimensional Hilbert spaces sup- 
porting irreps of SU(2) q labelled by j h with ji < (£ - 2)/2 (i = 1, 2, . . . , 2n). 
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• £?2n is the braid group on 2n strands whose generators g = {<7i, 02, ■ ■ ■ , &2n-i} 
and their inverses represent the input alphabet (we may add the identity element 
e € B-2n as null symbol). 

• Y = { accept, reject, e } is the output alphabet. 

• The transition matrices are expressed in terms of U(B2n)> denoting collectively 
the Kaul unitary representation matrices, while the projectors P(y) (y E Y) 
are defined as in the general case given in section 2.1. 

According to the above definitions, we provide the automaton with an input word w S 
Bm of length k (written in the alphabeth g by natural composition in Bin) 

w = < <\ ■ ■ ■ <?Z ; ff a ,es, ti = ±1 (18) 

and such that the (plat) closure w of the 2n-strand braid w gives the link L to be 
processed. Dropping for simplicity all the matrix indices, the unitary evolution of the 
automaton is achieved by applying the sequence 

U(w) = U(*H ) U(a^~_\ )■■■ U(al\ ) (19) 

to a start ket |0; 0) J in the odd-coupled basis ( fT~3T >. 

Whenever an odd-braiding a a = an+i (or (o^i+i) -1 ) occurs, the automaton one- 
step evolution upgrades the internal state by inserting the eigenvalue of the associated 
unitary U. On the other hand, when an even-braiding up = an (or (<72i) -1 ) must 
be implemented, the automaton has to change the parity of the internal state by means 
of a duality matrix (fT5T >. so that the effective transformation is given by the product 
U(ap) A[: :]. Since any duality transformation can be split into a sequence of basic 
duality matrices, we may look at the q-Qj symbol as representing an 'elementary', 
one-step evolution of the automaton. By resorting to standard results in graph theory 
it is possible to estimate how many q-Qj symbols are needed to decompose the most 
general q-recoupling coefficient (67), In the present case the upper bound can be 
expressed in terms of 2n, the braid index (or, equivalently, the number of strands of the 
input braid). On the basis of the above remarks, the time complexity function (num- 
ber of computational steps) for processing a braid-word of length n on the quantum 
recognizer A q is bounded from above by (2] 

k (N InN) where N = (2n - 1), (20) 

implying that the automaton processes efficiently such braids. 

Let us finally comment on the 'probability distribution' entering into the definition 
of a quantum automaton that recognizes a language in a probabilistic sense (end of 
section 2.1). On the basis of the expression of the colored link invariant given in ( TToT l 
and by comparison with the word probability of a quantum recognizer defined in <[8j>, it 
should be quite clear that the probability naturally associated with a link L processed on 
A q is the square modulus of its colored polynomial (note that the positivity conditions 
required in (O are always satisfied). 

In order to check this result in a concrete case, the explicit construction of the q— 
spin network automaton that recognizes the braid group language with a probability 
distribution given by the square modulus of the Jones polynomial is carried out in 
appendix B. 
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3 Efficient quantum algorithms for 3-manifold 
quantum invariants 

3.1 Colored framed links and 3-manifold quantum invariants 

The quantum invariants of 3-manifolds that we are going to discuss -within the math- 
ematical framework developed in 11301 . see also l68l - can be obtained as combinations 
of polynomial invariants of 'framed' unoriented links in the 3-sphere S 3 on the basis of 
theorem 1 stated below. It is worth noting that in the CS W environment the necessity of 
introducing framings is physically motivated by the requirement of general covariance 
of the quantized field theory (see e.g. ch. 3 of ll23lD . 

Loosely speaking, a framed oriented link [L; f] is obtained from a link L-thought 
of as made of knotted strings- by thickening its strings to get oriented 'ribbons' . If L 
has S knot components K\ , K2 , K$, for each K s we introduce another closed path 
K( oriented in the same way as K s and lying within an infinitesimal neighborhood of 
K s (knots and links are embedded in M. 3 or in S 3 ). The overall topology of the link 
is not modified, but for each K s we now have an extra variable t(K s ) telling us how 
many times the oriented ribbon is 'twisted'. Denoting by f = {/ s = n(s), n(s) E Z} 
(s = 1, 2, . . . , S) the framing of the link L, f s is the self-linking number of the band, 
or equivalently the linking number lk(K s , K[) = x(Ks, K[) between the knot K s 
and its framing curve K[ which winds n(s) times in the right-handed direction. The 
twist of the band t(K s ) is not independent from lk(K Sl K(), and the simplest choice 
we can made is to set 

t{K s ) = w{K s ), 

where w(K s ) is the writhe of the s-th component. Q 

The type of framing usually adopted is the 'vertical' framing, where the frame 
is thought to be placed vertically above the link diagram. Looking at Fig. |9), the 
framing can be represented by putting the writhe number w(K s ) = n(s) nearby the 
band or even by writing n(s) next to the string representing (a portion of) K s . 

The 3-manifolds we are going to consider are closed, connected and oriented, 
namely compact and without boundary. In general any such manifold Ai 3 can be pre- 
sented as the union of several components -endowed with 2-dimensional boundaries- 
glued together by suitable identification prescriptions on the points lying in their bound- 
aries. If the components are sewed together in a different way a topologically different 
manifold A4 3 may be obtained (recall that two manifolds are topologically equivalent 
iff they are homeomorphic, namely there exists a one-to-one continuous map between 

4 Given a link diagram D(L), namely a projection of the oriented link L onto a plane, we can define 
two numerical invariants associated with such diagrams. The writhe number w(Dl) is given by w(Dr / ) = 
e (p) where {p} are the crossing points £ D(L) and e(p) = ±1 according to whether there appears an 
over-crossing of the left strand over the right strand or an under-crossing (both strands are oriented upward). 
The linking number, defined for a link with more than one component knot, is defined for each pair of 
components (Ki,Kj) as lk(Kj, Kj) = w{Dj 1 ) — w{Dx i ) — w(Dk^ ). 

It can be shown that the writhe is a regular isotopy invariant for knots and link and the linking number is an 
ambient isotopy invariant for links. Recall that two links in R 3 (or S 3 ) are ambient isotopic if they can be 
continuously deformed one into the other. It can be shown that two links are ambient isotopic if and only if 
their diagrams are connected by a finite sequence of Reidemeister moves of type I, II, II. Regular isotopy is 
a restricted type of equivalence among links where the allowed Reidemeister moves are of type II and III. It 
is worth noting that the colored link polynomials )16t are invariants of regular isotopy while an associated 
ambient isotopy invariant can be obtained by multiplying J(L; j; q) by {g- 3 ™^)/ 4 /(g 1 / 2 - q' 1 / 2 )}, 
where w(L) is the writhe defined above. 
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vertical framing 



Figure 9: Graphical representations for the vertical framing. 

them). The new manifold A4 3 is said to be obtained from A4 3 by means of 'surgery'. 
The fundamental theorem that characterizes (equivalence classes of homeomorphic) 
3-manifolds reads [|69ll70l 

Theorem 1. Every closed, connected and orientable 3-manifold A4 3 can 
be obtained by surgery on an unoriented framed link in the 3-sphere S 3 . 

Following ll23l . let us illustrate in some details the (Dehn) surgery procedure referred 
to in the statement of the theorem in the simple case of a knot K (1-component link) 
C S 3 . The building blocks of the construction are solid tori, so let us begin the discus- 
sion by considering the 2-dimensional torus T 2 = S 1 x S 1 (the cartesian product of 
two unit circles). A point in T 2 can be specified by the coordinates [e 101 , e 102 ) in the 
complex plane and in particular a point lying on the longitude has coordinate (e 101 , 1) 
while a point in a meridian has coordinate (1, e 102 ). A solid torus V 3 is a 3-manifold 
homeomorphic to S 1 x T> 2 , where V 2 is the 2-dimensional disc. The boundary dV 3 
of V 3 is the torus T 2 . 

A framing on the solid turns is a particular homeomorphism h : S 1 x T> 2 — > V 3 . Given 
a framing h of V 3 , h(l x &D 2 ) is a meridian while h(S 1 x 1) is a longitude. 
A tubular neighborhood M of a knot K in S 3 is an embedding l : K x T) 2 —> S 3 such 
that t(x, 0) = x Vie <E K. The framing K$ of a knot K is defined as the framing of 
the tubular neighborhood of the knot. In particular, the standard framing K* of a knot 
K is such that the linking number x(K, K*) is equal to zero. 

Representing the unit disc T> 2 in the complex plane, the points of S 1 x T> 2 have coor- 
dinates (e 101 , re 102 ) with < r < 1. The self-homeomorphisms t of S 1 x T> 2 , defined 
explicitly by 

i±(e i8 \ re 102 ) = (e i9 \ re (e2±ei) ), (21) 
are the basic twist operations of S 1 x T> 2 . 

Given a tubular neighborhood J\f of a knot JC with standard framing / : S 1 x T> 2 — > A/", 
the right-handed and left-handed twist t± of Af are related to the t± by 

t± = fi±r 1 - (22) 

A De/m surgery performed along a (framed) knot K C S* 3 can be described as 
follows. 
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o 



1. remove the interior J\f of a tubular neighborhood Af of K (the resulting manifold 

o 

(S n A/") is the complement torus); 

o 

2. consider (5 \ A/") and J\f as distinct spaces; 

o 

3. glue back M and (S \ AT) by identifying the points in their boundaries through 

o 

a given homeomorphism h : cW — * d(S \ AO- 
The resulting manifold A^ is recovered by setting 



and it is completely specified by the knot K and by the chioce of the gluing homo- 
morphism h. Equivalently, the surgery is characterized by the knot K and by a closed 
curve 7 6 dAf representing where /i is the meridian of JV. 

Dehn surgery is a simple and constructive prescription which basically consists 
in removing and sewing back solid tori from the 3-sphere. However, since different 
surgery instructions may give rise to homeomorphic manifolds, it is crucial to define the 
equivalence relations that identify the surgery instructions providing the same (homeo- 
morphism class of) 3-manifold. Once these rules are taken into account, the classifica- 
tion problem for 3-manifolds can be actually reduced to the problem of classification 
of knots (links). 

The equivalence relations among surgery instructions yielding a same 3-manifold 
are topological operations on framed link diagrams known as Kirby moves. 
I move. The configuration described by the unknot U, with framing +1, enclosing 
n unlinked strands Ki, with framing fi, lying on a ribbon can be changed into the 
configuration where U is removed and the ribbon is twisted in the clockwise direction 
from below, see Fig. 1101 . The framing of the components Ki's must be changed 
according to 



II move. An unknotted link U with framing —1 can be removed without affecting the 
rest of the link, see Fig. fTTI ; 




(23) 



h 




(24) 




Figure 10: The first Kirby move. 
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Figure 11: The second Kirby move. 



III move. This is the 'inverse' of the first move and amounts to change the configura- 
tion containing U with framing —1 into the anti-clockwise twisted ribbon of AYs and 
without Here the framings change according to 

fi -> f'i = h + X^^U). (25) 

IV move. This is the 'inverse' of the second move. An unknotted link U with framing 
+1 can be removed without affecting the rest of the link. 

The extension of the results stated above in order to deal with 5J7(2) g -colored 
surgery links to be associated with 3-manifold invariants can be carried in a quite 
straightforward way. The construction developed in J2] |3] and summarized in section 
2.3 basically relies on Kaul unitary representation of (colored) oriented braids on In 
strands [28 1 and thus we need preliminarly to generalize such a representation to deal 
with unoriented braids and links. The goal can be achieved by suitably modifying 
the eigenvalues of the elementary braiding matrices in the odd-coupled basis |p; r) J 
introduced in dT3l according to ll29l 

U(a 2 i+i ) |p;r) (-'hi+i, 321+2,-) _ \ t | p; r ) ( '" 'hi+i, 321+2,- ) 
with A, = X(j 2l+1 , j 2l+2 ) = A,(j, f) = qHci+o,,-^ )/ 2) (26) 

where cj = j(j ' + 1) and cy = are the quadratic Casimir invariants associated 

with the irreps j, f . The ±1 in the exponent of the parameter q refer to left-handed 
(respectively, right-handed) half-twists in two parallel strands carrying the coloring 
j, j' and it can be easily checked that these eigenvalues do not depend on the orien- 
tations but only on the over/under-crossing features (the colored strands are ordered 
from left to right as happened for oriented braids, see dTOb and Fig. J3|). 
It is worth noting that the duality matrices (JT3J -needed whenever an even braiding 
U(a 2 i) has to be applied- are independent from orientations of the strands, so that they 
can be used in the present context with no further modification. 

According to Theorem 1, what we really have to handle are colored links in the 
vertical framing f and then the effect of adding or deleting a ±1 in the writhe of the 
link must be properly taken into account (this operation in the standard framing would 
not affect the topology of the link). Referring in particular to a j-colored framed 
unknot U, the associated link invariant turns out to be changed into 

J[U;j,+l;q] = q Cj J [U; j, 0; q] 

J[U;j,-l;q] = q~ Ci J [U; j, 0; q], (27) 



21 



where ±1 on the left-handed sides denote the ±1 vertical framing while in the right- 
handed sides there appears the knot invariant of the unknot in the standard 0-framing, 
whose numerical value is given by the ^-dimension [2j + l] q defined in (T% . The latter 
relations provide in practice the operatorial content of Kirby moves applied to framed 
colored links. Then the requirement of invariance under Kirby moves of the forthcom- 
ing 3-manifold topological invariants makes it necessary to compensate the effects of 
dZTt by resorting to properties of the so-called linking matrix. 

For a framed link [L; f] whose components K\ , K 2 , ■ • • , Kg have framings m , n 2 , ■ ■ ■ , 1 
respectively, the linking matrix is a symmetric matrix defined as 



X [L; f] 



/ m x(Ki,K 2 ) X (Ki,K 3 ) ■■■ x(K 1 ,K s )\ 

x(K 2 ,Ki) n 2 X (K2,K 3 ) ■■■ X {K 2 ,K S ) 



\ \;A-A:i ••• n s J 



(28) 



where \(Ki,Kj) is the linking number between the component knots Ki and Kj. 
The signature of the linking matrix, denoted by <r[L; f], is the difference between the 
number of positive and negative eigenvalues of X [L; f]. 

With these preliminary definitions, let us state the following theorem, the original proof 
of which can be found in |[l2"l[30l (see also [68 1). 

Theorem 2. For a closed, connected and oriented 3-manifold M\ ob- 
tained by surgery in the 3-sphere along an unoriented colored framed link 
[L; f , j] with S link components and for any fixed root of unity q = e' E + 5 
the quantity 

I [Ml ;f;q] = a- ff[L; fl ]T thx Hi* ■ ■ ■ J \L; f, j ; g] (29) 

{j} 

is a topological invariant of the 3-manifold endowed with the framing 
assignment f |f| 

3irik ,, . _ / 2 • tt(2j + 1) . _ 



Herea = ex Pi ^||y,^ = ^ ^ sin ] ,j = j 2 , . . . ,js) runover{0, |, . . . , 
and the summation is performed over all admissible colorings. J [L; f , j ; q] in the pre- 
vious expression is the 'unoriented' counterpart of the polynomial for the link L with 
coloring assignment j on its components and the summation is performed over all ad- 
missible colorings. 

It is worth noting that the presentation of the colored links used both in section 2.2 
(see (fTol l) and in ( f26b is slightly different from the presentation given in Theorem 2 
above, even though we keep on using the same notation j for the colorings. In the latter 
case a coloring is assigned to each of the S link components, while in the former we 
label the 2n strands of the associated braid with n colors (see e.g. Fig. 01). 



5 This means in practice that, on applying Kirby moves I-IV to [L; f] C S 3 , the value of the invariant 
does not change, namely depends only on the homeomorphism class of the 3-manifold. Note however that 
we are not in the presence of a complete 3-manifold invariant since there exist topologically distinct mani- 
folds with the same I. 

The extension of the theorem to deal with surgery operations performed on manifolds topologically different 
from S 13 and to situations in which g-deformations of other semisimple Lie groups are involved can be found 
in the reference quoted above. All such invariants are collectively referred to as 3-manifold 'quantum invari- 
ants' or even as 'state sum models' by noticing that an expression like {29} can be interpreted as a partition 
function over suitably weighted 'states' represented by the link polynomials with different colorings. 
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However, what really matters is the fact that both presentations of the link L give rise 
to the same invariant, as could be proved by exploiting the properties of the SU (2) q — 
representation ring. From the computational viewpoint this twofold choice does not 
matter as well because we could efficiently implement the transformation from a given 
link diagram to any associated closed braid on a classical computer (see e.g. sec- 
tion 2 of lETTl for a discussion on classical algorithmic questions about braids and 
links). Expression d29l makes it manifest the overall dependence on the framing in 
the factor a^ a ^ L ' f l, but in what follows we are going to express the link polynomials 
J [L; f , j ; q] as expectation values of Kaul unitary representation [29 1 worked out for 
an unoriented 2n-strands braid the plat closure of which gives the link under examina- 
tion. 

3.2 Quantum algorithm for approximating 3-manifold invariants 

In this section we shall describe the quantum algorithm for computing the colored 
polynomial J [L; f,j; q) and the associate invariant X \M\ ; f; q] defined in j29l. This 
algorithm is an extension of the quantum algorithm proposed in [UJEI which efficiently 
approximated the value of the colored Jones polynomials ([Tol l. 

As anticipated in section 1.2 we need in the present context the notion of addi- 
tive approximation introduced [14] (see also [16]). Given a normalized function g(x), 
where x denotes an instance of the problem, we have an additive approximation of its 
value for each x if we can associate with g{x) a random variable Z such that 

pr { \g(x) - Z | < 77} > 3/4 , (30) 

for any 77 > 0. Moreover, the time needed to achieve the approximation must be 
polynomial in the size of the problem and in r\~ x . Then the problem we are interested 
in can be stated formally as follows. 

Approximating 3-manifold invariants. 

Given a framed link [L; f] with component knots K\, K2, ■ ■ ■ , K$, fram- 
ing f = (rti, ri2, . . • , lis), a positive integer k, and a set of allowed col- 
ors j = {0, 1/2, . . . , fc/2}, we want to sample out a random variable Z 
-representing an additive approximation of the value of the normalized 3- 
manifold invariant X \M\ ; f ; q] evaluated at q = e 7 ^ - in such a way 
that the following condition holds true 

pt{\l[M 3 L ;f;q] - Z \ < 77} > 3/4 . (31) 

The size of the problem is expressed in terms of the number of crossings k 
of the surgery link L and by the number of strands of the associate braid, 
(as it happened for colored link polynomials), but we shall need to handle 
properly the whole set of allowed colorings j and the framing f as well@ 

As anticipated in section 2.1, the quantum algorithms for evaluating (additive ap- 
proximations of) topological invariants in the framework of quantum CSW theory are 
based on a two-level procedure outlined already in [3 1 in connection with colored link 
polynomials and improved here for 3-manifold invariants. The rationale underlying 
our procedure is briefly summarized below while technical details are developed in the 
rest of the section. 

6 From now on we agree that the invariant is normalized by the product of the g-dimensions associated 
with the link components, namely by the factor nf=i Pii + 1]<j- 
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(A) Within the computational model of the q-deformed spin network recognizer A q 
(section 2) both topological and field-theoretic data -encoded into a framed col- 
ored link [L; f, j ] associated with a unitary braiding operator in the Kaul rep- 
resentation [29 1- are efficiently processed. The estimate of the time complexity 
function (number of computational steps) required to complete the calculation 
given in ( f2Qb still holds true. 

(B) By resorting to standard quantum circuit model techniques and related approxi- 
mation schemes it can be shown that 

a) the start state \ip) x of the recognizer A q needed for processing the 3- 
manifold invariant I [M\ ; f ; q] is efficiently encoded into a qubit register; 

b) the braiding operator associated with the framed link [L; f , j ] -already split 
into a sequence of 'elementary' braiding and duality transformations on the 
basis of the recognizer design- can be efficiently compiled on a standard 
quantum circuit applied to the start qubit register; 

c) by resorting to the Hadamard test [15] -an efficient sampling procedure 
which provides the expectation value of an unitary operator on a (qubit) 
state- it is possible to estimate the value of the invariant X \SA\ ; f ; q] from 
a series of measurements on an ancilla qubit coupled to the start state. 

This two-level computational process for the approximation of 3-manifold in- 
variants in the sense of OTJ is thus efficient with respect to the standard model 
of quantum computation and not simply from the viewpoint of the g-deformed 
automaton model. 

(C) The expectation value of the braiding operator evaluated on the start state and 
sampled as described above, can in turn be related to a suitable probability dis- 
tribution on the language(s) recognized by A q . In view of the properties of such 
distributions (cfr. the concluding remarks of sections 2.1 and 2.2 and appendix 
B), the whole procedure can be reinterpreted by saying that the automaton rec- 
ognizes the language of the braid group with a probability distribution given by 
the square modulus of the (normalized) invariant 1 [M,\ ; f ; q]. 

Given an unoriented link L presented as the plat closure of a 2n-strands braid with 
a fixed set of colorings j = ji, j'2, . . . , j2n ( see section 2.2 for the oriented case, in 
particular Fig. J4)), let us denote again the conformal block odd-coupled basis of the 
boundary WZW theory by |p; r ) J as in (TT~3T > and Fig. J5). Here p = po, . . . , p n -i and 
r = ro, ... , r„_3 and both j and the intermediate quantum numbers p, r take values 
in the collection {0, 1/2, ... , fc/2}, bounded form above by the coupling constant k of 
quantum CSW theory. 

Each spin quantum number can be encoded into a qubit register made of [ log 2 (k + 1 ) ] 
qubits, where \ x~\ denotes the smallest integer > x, see Fig. 0121 . 

Since an element of the basis is specified by (An — 3) quantum numbers (j,p and 
r), we need 

(4n-3) x pog 2 (fc + l)l (32) 

qubits to encode one basis vector. The ordering on such quantum register is shown in 
Fig. 1131 . where the three sets of quantum numbers are associated with a j-register, a 
p-register and a r-register respectively. 
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#flog(k+1)l 

Figure 12: The i-th quantum number is encoded into a register made of [ log (fc + 1) ] 
qubits. 

00 -000- 000 -g 

j-register p-register r-register 

Figure 1 3 : The register of the quantum circuit encoding an odd-basis state | p ; r ) J . 

The unitary braiding operator associated with the 2n-strands unoriented braid un- 
der consideration can be decomposed -following the scheme explained in section 2.2 
and updating the representation according to |29|- into 'elementary' odd braidings and 
duality trasformation. 

The elementary braiding matrices are diagonal in the odd-coupled basis, so that their 
action can be easily implemented on the quantum register of Fig. lfT3l . The quantum 
gate realization of U(a2i+i) is simply the identity matrix on the (p, reregisters, while 
the j-register is modified by the action of a SWAP gate and a phase gate with a phase 
factor given by the eigenvalue A(j 2 /+i , 321+2) in equation d26ill 7 l 

The duality transformation introduced formally in dl~5b must be applied whenever 
an even elementary braiding matrix J7(ct2/) is encountered in order to recover the cur- 
rent, odd-coupled basis. The associated (g-deformed) 3nj coefficient can be split into 
a sequence of 'elementary' duality transformations according to 




i(q; s ) 
L (p;r) 



3\ 32 
33 k 

32n-\ 32n 
n-2 



£l...t n _2 2—1 

n-2 

W_ Ariqi+i 



n-i 32i+\ 
]2i+2 n 



A 



tiSi — i 



1=0 



tl 321+2 
321+3 U+l 



'Recall that a SWAP acting on two qubits \x), \y) is the operation 

SWAP: \x)\y)^ \y)\x), 

which corresponds to the matrix 

/ 1 \ 
10 
10 

y 1 ) 



SWAP : 



U-l Qi 
Si 32r, 



(33) 
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Figure 14: Circuit diagram of the unitary gate associated to an odd braid generator. 



(see Fig. Q of section 2.2 for the graphical representation of such decomposition). 
Each symbol in the latter expression is the matrix form of a q-6j coefficient (see Fig. 
ID), namely 



AJ12 



Jl 32 
33 3 



= (-)Wi+*+i.-hO ([ jy la + 1] [2m + 1] Jl * J12 

1 J3 J 323 



1 

(34) 

where the labels of the quantum numbers has been slightly changed to comply with the 
following standard explicit expression of the q-6j 



Jl 32 3X2 
33 3 323 

]2, ii2)A(j 3 , j, j 12 ) A(ji, j, j 23 ) A(j 2 , 33, 323) 

(-)'[* + 



E 

z>0 



[z - 3i - 32 - ji2]ql[z - 33 -j- 3n]q\[z - 3i -3- 33W- 



[z - 32 - 33 - 323]q ] -\jl + 32 + 33 + 3 ~ z ] q \ 



[jl + 33 + 312 + 323 - z] q \[j 2 +j+ 312 + 323 ~ z] q \ 



(35) 



The symbols A(. . . ) are combinatorial factors involving the factorials of q-dimensions 
(see appendix I of (28]), where the factorial of a g-number [x] q is defined as [x] q \ = 
[a;] g [s-l] 9 ...[2] 9 [l],. 

For each choice of the entries (in the allowed set {0, 1/2, . . . , fc/2}) the latter power 
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series is actually a summation over the finite set of (integer and half-integers) z's that 
yields non-negative quantum integers. 

From the remarks above it should be clear that the problem of efficiently imple- 
menting d33l ). namely the most general change of basis, is equivalent to the simpler 
problem of efficiently compiling a sequence of q-6j symbols or elementary duality 
matrices d34l i. But any such coefficient (i.e. the symbol with a fixed set of entries) can 
be easily and efficiently evaluated with a classical computer owing to the finiteness of 
its explicit expression 1351 . 

On the other hand, it is necessary to make explicit the action of the matrix d34l) on the 
qubit register, namely on an Hilbert space of dimension 2 F log ( fc +t)l for each admissible 
set of entries. @ 
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Figure 15: The multiplexor gate representing the action of a q-6j transformation. 



The circuit realization is shown in Fig. 1151 . where the j-register acts as a control 
register on the qubits of the p(r)-registers involved in the transformation. The trans- 
formation A can then be thought of as a gate with a 'block structure', or a multiplexor 
associated with the block structure depicted in Fig. 11161 . Each block corresponds to 
a particular configuration of the j-qubits, and the matrix element inside the block are 
q-6j, up to suitable factors (equation (|34|>). 

On the basis of the decomposition in (l33l > we realize that the allowed elementary 
duality matrices are always parametrized by the set j of those quantum numbers which 
remain unchanged when they are applied to the proper registers. The crucial remark 
here consists in noticing that the dimension of these matrices is independent of the 
size of the problem, determined by the index of the braid group and the number of 
crossings. Since there exist efficient methods to approximate unitary matrices of a 
given dimension [71|, a sequence of universal gates can be always worked out that 
efficiently approximate every q-6j as well (see 11721 for more details). 

The number of elementary duality transformations needed to decompose a general 
duality transformation is (2n — 3), linear in the size of the problem under consid- 
eration. The action of a q-3nj recoupling transformation on the register of Fig. fT3l 
is shown in Fig. ifTTl . Note however that, once all the (2n — 3) gates represented by 
q-6j have been applied, it is necessary to swap some of the qubits in order to recover 
the proper order in the register, see Fig. lfT~8l . 

Upon applying a q-3nj recoupling transformation we end up in the even-coupled 
basis |q;s)-> (see (fT4l i and Fig. J6)) which diagonalizes the even braiding matrices 
U(a2i) of the Kaul representation. Their action on the (even) qubit register can then be 

8 The 'classical' 6j'-symbol can be expressed in terms of a series similar to (35), but the encoding of such 
symbols on qubit registers would depend explicitly on the values of the entries, so that the computation is 
not obvioulsy efficient with respect to the size of the input data, see footnote 2 in section 2.1. 
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Figure 16: The block-decomposition of the matrix corresponding to a q-6j transfor- 
mation. 
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Figure 17: The circuit realization of a general q-3nj recoupling transformation. 
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Figure 18: The sequence of SWAP gates needed to reordering the qubit register. 
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implemented by paralleling the procedure described in the odd case. 

The further step consists in going back to the |p;r) J basis, which can be achieved by 

means of the finite sequence of q-6j transformations shown in Fig. (19), eventually 

completed by a suitable sequence of SWAPS to recover the initial ordering of the start 

register. 

11) a a 
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|5) 1 1 « 
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I m 1 .11 I 111 TT^ 

I 10) 

11) 

|2) 
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|n-1) 

|1) 



,M 

Figure 19: The inverse of the circuit of Fig. H17I , 

This completes the analysis of the quantum circuit which implements efficiently 
the Kaul representation 1 29 1 on the Hilbert space spanned by the odd-coupled basis of 
conformal blocks (point b) at the beginning of this section). 

For what concerns the preparation of the start state according to point a), notice 
once more that we have considered so far the case of an unoriented link [L; j] (plat 
closure of a 2rt-strands braid) with a fixed coloring set j. As already pointed out in 
the case of oriented links (see section 2.2, (TToT i and appendix B) the start state of the 
automaton calculation is a singlet vector in the odd-coupled basis 

[L;j] <-> |0;0) J ; p = p , . . . ,p n _i = , r = r , . . . , r„_ 3 = 0, (36) 

which can be efficiently encoded into the qubit register of Fig. |[T3"1 as a particular 
case of the construction carried out above. Then the polynomial for an unoriented, 
framed link J [L; f, j ; q] can be recovered as the (trace of the) expectation value of 
the associated braiding operator in the updated Kaul representation. However, the 
start state \ijj) 1 of the recognizer A q needed for processing the 3-manifold invari- 
ant I ; f ; q] must be prepared in order to comply with the summation over the 
colorings in ( f29l . The goal can be achieved in a quite staightforward way by resorting 
to quantum parallelism, as described in the following. 

As before, we first split the qubit register into three ordered parts, the j-, p- and 
r-registers. Then we initialize the j-register into the following weighted superposition 
of states 

l^) J = E |J>, (37) 
{j} 

with /ij = (ttjj • • • fij a (S in the number of link components and the coefficients fi's are 
defined as in Theorem 2 of section 3.1). Since the link L can be considered as the plat 
closure of a braid word w, the quantum circuit needed to encode the start state can be 
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implemented by the following composition of unitary gates 

S J (w) U mP ® R (w) : J <g> |0; 0) h-> |^) J ® |p'; r') P ® P (38) 

where S J (w) is the SWAP gate on the j-register and |p'; r') is the state obtained from 
the application of the circuit U(w). 

The following sequence of equalities shows that this quantum circuit actually gives 
the required 3-manifold invariant 

p ^(0;0| ® J (iP\S J (w)U J ® p ® R (w)\ij) J |0;0) p ® p = 
= P ® R (0;0\ ® J (^| • |^) J ® U p ® R {w)\0;0) p ® R = 
= £ W <0;0|tf p ®>)|0;0> = 

{j} 

= I[M|;f; 9 ]a^ f l. (39) 

Recall that the signature of the linking matrix <r[L; f] can be easily computed classi- 
cally once the linking matrix (f28b (which is part of the topological input data of the 
algorithmic problem) is given. 

Finally, the efficient sampling of the value of the invariant referred to in point c) 
can be carried out by resorting to the techniques of lfT5l[l6l . as already done in J51 for 
the case of colored link polynomials. 

4 Concluding remarks 

We have shown that all the significant quantities -partition functions and observables- 
in SU(2) quantum CSW theory can be efficiently approximated by quantum algo- 
rithms at finite values of the coupling constant k, reflecting the intrinsic field-theoretic 
solvability of this theory. The efficiency of the quantum algorithms is proved on the 
basis of a two-level computational scheme which relies on the implementation of uni- 
tary representations of the braid group proposed in |28, 29]. In particular, as shown in 
section 3.2, the representation needed to handle 3-manifold quantum invariants can be 
efficiently compiled on a quantum circuit equipped with a suitable start qubit register, 
namely within the 'standard' model of quantum computation. 

The relevance of our result in connection with the algorithmic complexity of quantum 
CSW theory has been extensively addressed in section 1 .2, while further developments 
and applications to geometry and quantum gravity models had already been discussed 
in 1.3. 

More generally, it would be interesting to improve our approach, on the one hand, 
by handling quantum topological invariants arising from CSW theories with arbitrary 
semi-simple Lie groups ll73l and, on the other, by exploring the quantum computa- 
tional complexity of solvable models in statistical mechanics (74). 

Finally, let us comment in some more details the issue concerning the model of 
quantum computation adopted . As pointed out in (A) and (C) of section 3.2, a central 
role in our procedure is played by the quantum recognizer A q able to process efficiently 
the language generate by the braid group with transition matrices given by Kaul unitary 
representations (section 2.2 and 3.2) and probability distributions associated with quan- 
tum topological invariants. It may be tempting to proceed without this step, processing 
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directly the unitary representations within the quantum circuit scheme of computation. 
However, we proved in section 3.2 that the basic morphisms of the q-tensor category 
( D\(SU (2)«j ) ; 1Z ; T ) on which the recognizer is modelled can be efficiently com- 
piled and approximated on a quantum circuit (in particular the implementation of a 
q-Qj transformation is independent of both the input size of the algorithmic problem 
and on the values of its entries; see also footnotes 2, section 2.1 and 8, section 3.2). 

This achievement, quite remarkable as it is by itself, opens as well the further possi- 
bility of looking at the q-spin network simulator as the fundamental model of computa- 
tion for a wide range of algorithmic problems in geometric topology and group theory. 
According to the quantum recognizer definition given in section 2.1, the specific prob- 
lems that can be dealt with it require the selection of suitable start states -that must be 
efficiently encoded into qubit registers- and of (possibly constrained) sets of transition 
matrices. We conjecture that the approximate evaluation of invariants for (colored) tri- 
angulations in dimension 2 and 3 should require minor modifications of the scheme 
employed in this paper for processing quantum invariants of links and 3-manifolds. 
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Appendix A. The spin network quantum simulator 



The universal model of quantum computation proposed in ATI l42l l43l relies on the 
(re)coupling theory of SU(2) angular momenta l56l l60l . It can be thought of as 
a generalization of the standard quantum circuit model in which the computational 
Hilbert spaces are binary-coupled eigestates of N = (n + 1) SU (2)-angular momenta 
(whose quantum numbers range over {0, 1/2, 1, 3/2, . . .}) and unitary transformations 
('gates') are expressed in terms of recouplingcoefficients (3nj symbols) connecting 
pairs of inequivalent binary coupling schemesO 

The architecture of the spin network is modelled as an SU (2) fiber space structure 
over a discrete base space V 

(V, C 2J+1 )„ (40) 

which encodes all possible computational Hilbert spaces as well as all gates for any 
fixed number N = n + 1 of incoming angular momenta (see appendix A of [42 1 
and [43 1 for more details). The base space V = {^(b)} represents the vertex set 
of a regular, 3-valent graph, the so-called twist-rotation graph [67 1 &„(V,E) with 
cardinality \V\ — [2n)\/n\, i.e. the quadruple factorial number. E, the edge set of the 
graph, will be associated with permitted transformations between pairs of verices as 
described below. 

There exists a 1 : 1 correspondence 

Wb)}^R 7 (f>)} (41) 

between the vertices of (& n (V,E) and the computational Hilbert spaces of the simu- 
lator, where b denotes the binary bracketing structure that we are going to describe. 
For a given value of n, H-n(b) is the simultaneous eigenspace of the squares of In + 1 
Hermitean, mutually commuting angular momentum operators 

Jl) J2, J3) • • • j Jn+l = {J?:}; Jl + J2 + J3 + •■■ + Jn+1 = J 

and K x , K 2 , K 3) . . . , K„-i = {K,J (42) 

together with the operator J z (the projection of the total angular momentum J along the 
quantization axis). The associated quantum numbers are ji, fa, . . . , j n +i ; J; k\, k% , . . . , 
k n -i and M, where — J < M < in integer steps. If W 1 (8 W 2 <g> • • • ®W n ® W n+1 
denotes the factorized Hilbert space, namely the (n + l)-fold tensor product of the 
individual eigenspaces of the (J,) 2 's, the operators K^'s represent intermediate angu- 
lar momenta generated, through Clebsch-Gordan series, whenever a pair of Jj's are 
coupled. 

As an example, by coupling sequentially the Jj's according to the scheme (• • • ((Ji + 

J2) + J3) H h Jn+l) = J - which generates (Ji + J 2 ) = Ki, (Ki + J 3 ) = K 2 , 

and so on - we would get a binary bracketing structure of the type (•••(( (W 1 ® W 2 )ki 
®W 3 )k 2 ® ' " " ® T~L : '" +1 )k rl - 1 )j, where we add an overall bracket labeled by the quan- 
tum number of the total angular momentum J. Note that, as far as jj's quantum num- 
bers are involved, any value belonging to {0, 1/2, 1, 3/2, . . .} is allowed, while the 
ranges of the fc^'s are suitably constrained by Clebsch-Gordan decompositions (e.g., 
if (Ji + J 2 ) = Ki => \ji — J2I < fci < ji + 32)- We denote a binary coupled basis 

'The model can be extended to include Wigner rotations in the eigenspace of the total angular momentum, 
see section 3.2 of 1421 . 
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of (n + 1) angular momenta in the JM-representation and the corresponding Hilbert 
space appearing in ( |4TT > as 

{ I [ji, 32, is, ■ • ■ , jn+if ; fe? , k b 2 , . . . , k b n _ t ■ JM), ~J < M < J} 

= Hi (fa) = span { | fa ; JM) n } , (43) 

where the string inside j%, ja, . . . , j n +i] b is not necessarily ordered, fa is the short- 
hand notation for the current binary bracketing structure and the fe/j's are uniquely as- 
sociated with the chain of pairwise couplings selected by fa. 

For a given value of J each H„(E>) has dimension (2 J + 1) over C and thus there exists 
one isomorphism 

H J n {b) ^ b C 2J+1 (44) 

for each admissible binary coupling scheme b of (n + 1) incoming spins. It is worth 
stressing that such isomorphic spaces are physically inequivalent as far as they are 
associated with different schemes of (binary) interactions. 

The vector space C 2J+1 is interpreted as the typical fiber attached to each vertex v(b) S 
V of the fiber space structure d40b through the isomorphism (f44b . 

For what concerns unitary operations acting on the computational Hilbert spaces 
d43l . it can be shown [56 1 that any 3nj symbols of SU(2) can be splitted into 'el- 
ementary gates' represented by Racah and phase transformations (in the categorical 
language of section 2.1 they are referred to as the basic morphisms introduced in (01 
and (O, respectively). A Racah transform applied to a basis vector of the type d43l is 
defined formally as 

F : \...{(ab) d c) f ...;JM) ^ | ... ( a (bc) e )/...; JM), (45) 

where we are using here Latin letters a, b, c, . . . to denote both incoming (ji 's in the 
previous notation) and intermediate (fc/j 's) spin quantum numbers. The explicit ex- 
pression of d45l l reads 

\(a(bc) e ) f ;M) 

= ^(_l)^+c +/[(2d+1)(2e+1)] i/ 2 | a b d | K{ab)dc)f . Mh (46) 

where there appears the Racah- Wigner 6j symbol of SU (2) and / here plays the role of 
the total angular momentum quantum number. Owing to the Wigner-Eckart theorem, 
the magnetic quantum number is not affected by such transformation, and the same 
holds true for a general, 3nj recoupling coefficient. Recall also that the square of 
the 6j symbol in (|46| | represents the probability that a system prepared in the state 
| ( (a b)d c) f ; M) will be measured in the state | (a (b c) e )/ ; M). 
A phase transform on a basis vector (l43l is defined as 

R : \...(ab) c ...;JM) ^ \ . . . (ba) c . . . ; JM), (47) 

and esplicitly reads 

\...{ab) c ...; JM) = (-)-+ b - \ . . . (ba) c . . .; JM), (48) 

The edge set E — {e} of the twist-rotation graph <8 n (V, E) is a subset of the 
Cartesian product (V x V) where an (undirected) arc between two vertices v(b) and 
v(b') 

e(fa,fa') = (v(b),v(b')) e (VxV) (49) 
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exists iff the underlying Hilbert spaces are related to each other by an elementary uni- 
tary operation of the type d45b or ( |47T >. Note also that elements in E can be considered 
as mappings 

(V x C 2J+1 )„ — > (V x C 2J+1 ) n 

( V (b),ni(b)) -> (v(b'),ni(b')) (50) 

connecting each given decorated vertex to one of its nearest vertices and thus define 
a 'transport prescription in the horizontal sections' belonging to the total space (V x 
C 2J+1 )„ of the fiber bundle gOj. 

The structure of the graph <8 n (V, E) in the case of (n + 1) = 4 incoming spin variables 
a, 6, c, d is shown in Fig. 1201 (such a combinatorial pattern encoding both binary- 
coupled Hilbert spaces and trasformations among them was used for the first time in 

E3). 

a c d b 

.it. ,t„ i 

\y , w v V 



a d b c a d b c abed abed c da b c da b 



dead * ' « « ^ 



3 c o Z. <i 



\ 



b dc a b dc a 



Figure 20: A portion of the spin network graph 183 (V, -E). The vertices are labelled 
by rooted binary trees encoding the combinatorics of the binary coupled computational 
Hilbert spaces while the edges represent Racah transformations and the dashed ones 
are phases. 



The crucial feature that characterizes the graph <S n (V, E) arises from the compati- 
bility conditions relating the basic morphisms d45l l and ( f4Tb (referred to as the exagon 
and pentagon relations, see section 2.1). The Racah identity and the Biedenharn-Elliott 
identity together with the orthogonality conditions for 6j symbols (see e.g. [60 1 for 
their explicit expressions) ensure that any simple path in ® n (V,E) with fixed end- 
points can be freely deformed into any other, providing identical quantum transition 
amplitudes at the kinematical level. 

For what concerns the q-deformed spin network modelled on the g-tensor category 
( 9K(SU (2) q ) ; 1Z ; T ) defined in (O of section 2.1, we omit here all techical details 
and refer the reader to section 4 of [ 3 1 (see also the appendix of [ 28 1 for both definitions 
and notations). As already pointed out, the basic morphism T is implemented in this 
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case by means of a q-6j symbol defined in rt35l l while the braiding 1Z is to be associated 
with the over- or under-crossings of two contiguous strands belonging to (colored) 
braids (see e.g. the eigenvalue equation (l26l l of the odd braiding operators employed in 
section 3.2). 

A pictorial representation of ^-deformed spin network graph is shown in Fig. ETI . 




Figure 2 1 : A portion of the q-deformed twist-rotation graph (©3 (V, E)) q : with respect 
to the previous figure, each phase trasformation has been splitted in order to make 
manifest the non trivial braiding features. 



Appendix B. Automaton calculation 

The Jones polynomial ifTTI is a particular instance of colored link polynomial where the 
labels of the conformal block basis correspond to the fundamental, ^-irrep of SU (2) q . 
This allows us to simplify the notation for the states of the recognizer A q defined in 
section 2.2 by setting 

|p;r) J -> |p;r), 

since the coloring assignment {j} is always a string of 2n |— spins. The start state of the 
automaton is chosen to be |0; 0), i.e. all the internal labels are equal to 0, as depicted 
in Fig. 1221 . 




00 



Figure 22: The start state of the g-spin network automaton with probability distribution 
associated with the Jones polynomial. 

As discussed at the end of section 2.2, let us provide the automaton with an input 
word w of length k corresponding to an unitary evolution U (w) expressed in terms of 
the sequence given in jl9[ . 



35 



In order to complete the definition of A q we need explicit expressions for P (ac- 
cept) and P (reject) (see section 2.1). We choose the following 



P(accept) = |0;0)(0;0| , (51) 

P(reject) = I - |0; 0)(0; 0| . 

The 'a priori' probability distribution for the language generated by the braid group 
that we choose is the square modulus of the Jones polynomial of the plat closure w of 
the braid w, namely 

PrH = \J(w; q)\ 2 = V(w; q) , (52) 

where q is the root of unity at which the polynomial is evaluated. Using the properties 
of the Kaul representation (see (fTol l) and (Bit there follows that 

|Pr(iw)- (0; 0]^ (w)P(accept)U(w)\0;0) \ = 

\Pr(w) - (0; 0\U*(w)\0; 0)(0; 0|J7(«;)|0; 0)| = 
Pr(io) - |(0;0|C/(u;)|0;0)| 2 

V(w,q)-\(0;0\U(w)\0;0)f =0. 

Thus we have shown that the spin network quantum automaton recognizes 'exactly' 
(namely with a word-probability treshold 6 — 0, see (O in section 2.1) the braid group 
language according to the 'Jones probability distribution'. 

A similar result holds true for the families of automata A q parametrized by (ji , fa , 
■ ■ ■ jj2n) mat recognize the braid group language £?2n with a probability distribution 
given by the colored polynomial ( [Tol l. 
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